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We have calculated the large-q series of the energy cumulants, the magnetization 
cumulants and the correlation length at the first order phase transition point both 
in the ordered and disordered phases for the g-state Potts model in two dimensions. 
O^ , The series enables us to estimate the numerical values of the quantities more pre- 

cisely by a factor of lO'^ — 10* than the Monte Carlo simulations. From the large-q 
series of the eigenvalues of the transfer matrix, we also find that the excited states 
form a continuum spectrum and there is no particle state at the first order phase 
^ ' transition point. 

vn ■ 

f^ ■ 1 Introduction 

t^ ! 

f— ^ I In many of the physical systems that exhibit the first order phase transition, the 

f^ . order of the transition changes to the second order by changing the parameter 

^D \ of the system. It is important to know how the quantities at the first order 

phase transition point diverge when the parameter approaches the. point at 

which the order of the transition changes. The q-state Potts modetl'B in two 

dimensions gives a good place to investigate this subject. It exhibits the first 

flj ' order phase transition for g > 4 and the second order one for (? < 4. Many 

<~i . quantities of this model are knowELCxactly at the phase transition point (3 = f3t 

for q >.|4|^including the latent heau and the correlation length in the disordered 

phasepuo which vanishes and diverges, respectively, in the limit q —> A+. One 

^\ ' the other hand other physically important quantities such as the specific heat, 

S . the magnetic susceptibility, and the correlation length (in the ordered phase) 

at the transition point, which also diverge as g — > 4+, are not solved exactly. 

Here we calculate the large-g expansion series of the energy cumulants in- 
cluding the specific heat and the magnetization cumulants including the mag- 
netic susceptibility in both the phases and the correlation IcngjLh in the ordered 
phase at the transition point using the finite lattice methodQa'El Obtained long 
series for the energy and magnetization cumulants give the estimates of the 
quantities that are more precise by a factor of 10^ — 10"' than the Monte Carlo 
simulations. Especially its estimates are within an accuracy of 0.1% at g = 5, 
where the correlation length is as large as a few thousands of the lattice spac- 
ing. Bhattacharya et alcB made a stimulating conjecture on the asymptotic 
behavior of the energy cumulants at the first order transition point; the relation 



between the energy cumulants and the correlation length in their asymptotic 
behavior at the first order transition point for g — > 4+ will be the same as their 
relation in the second order phase transition for q = A and (3 —> Pt, which is 
well known from their critical exponents. The obtained series enables us to 
confirm the correctness of the conjecture. 

As for the correlation length at the first order phase transition point, the 
resultsjsf the Monte Carlo simulatiorO and the density matrix renormalization 
groujtJ indicate that the correlation lengths are very close to each other in 
the ordered and disordered phases for q > 10. On the other hand, at the 
second order phase transition point(<7 < 4) their ratio is known to be 1/2. 
It is interesting whether the ratio is exactly equal to unity, remains close to 
unity, or approaches 1/2 when g ^ 4_|_. To investigate it we calculate the 
first few terms of the large-g expansion for the eigenvalues of the transfer 
matrix and find that from the second largest to the 7V-th largest eigenvalues 
with N the one-dimensional size of the lattice make a continuum spectrum 
in the thermodynamic limit both in the ordered and disordered phases. We 
also calculate the long series of the second moment correlation length in both 
the phases, which serve to investigate the behavior of the spectrum of the 
eigenvalues of the transfer matrix in the region of q close to 4. 



2 Finite lattice method 

Here we use the finite lattice methodPISfl to generate the large-q series for 
the Potts modpL instead of the standard diagrammatic method used by Bhat- 
tacharya et afll3 The finite lattice method can in general give longer series 
than those generated by the diagrammatic method especially in lower spatial 
dimensions. In the diagrammatic method, one has to list up all the relevant 
diagrams and count the number they appear. In the finite lattice method 
we can skip this job and reduce the main work to the calculation of the ex- 
pansion of the partition function for a series of finite size,Iattices, which can 
be done using the straightforward site-by-site integratioioilj without the dia- 
grammatic technique. This method has been used mainly to generate the low- 
and high-temperature series in statistical systems and the strong coupling se- 
ries in lattice gauge theory. We note that this method is applicable to the series 
expansion with respect to any parameter other th an , t he temperature or the 
coupling constant. Using this method we calculateai30 the series for the n-th 
energy cumulants (n = — 6) to z'^^, n-th magnetization cumulants {n = 1 — 3) 
to z'^^ , and second moment correlation length to z^^ with z = l/i/g. 



Figure 1: Fade approximants of _Fi ' £'' at g = 4 plotted versus p 




3 Energy cumulants 

The latent heat C at the transition point are known to vanish at g ^ 4+ as 

with X = exp [tt"^ /29) and 2 cosh 6* — ^. Bhattacharya et al.'s conjecture says 

(n) 

that the n-th energy cumulants F^ J at the first order transition point (3 = f3t 
will diverge for g ^ 4+ as 



FP,{-1)"F^'''') r^aB"' 



r(|] 



3n/2-2 



(1) 



The constants a and B in Eq. ffl) should be common to the ordered and disor- 
dered phases from the duality relation for each n-th cumulants. 

If this conjecture is true, the product ^(")£3ri-4 ^g ^ gniooth function of 9, 
so we can expect that the Pade approximation of F^'^'C^ will give convergent 
result at p = 3n — 4. It has been examined for the large-g series obtained by 
the finite lattice methods for n = 2, • • • , 6 both in the ordered and disordered 
phases, which in fact give quite convergent Pade approximants for p — in — A 
and as p leaves from this value the convergence of the approximants becomes 
bad rapidly. An example can be seen in Fig. 1 for n=2 in the disordered phase. 
We give in Table 1 the values of the specific heat C — PfF'-^^ evaluated from 
these Pade approximants for some values of q. These estimates are three or 
four orders of magnitude more precise than (and consistent with) the previous 
result for q > 7 from the large-g expansion to order z^° by Bhattacharya et alt3 
and the result of the Monte Carlo simulations for (7 == 10, 15, 20 carefully done 
by Janke and KappleK3 as in Table 2. What should be emphasized is that we 
obtained the values of the specific heat in the accuracy of about 0.1 percent at 



Table 1: The specific heat for some values of q. The exact correlation length is also listed. 



q 


Cd 


Co 


^d (exact) 


5 


2889(2) 


2886(3) 


2512.2 


6 


205.93(3) 


205.78(3) 


158.9 


7 


68.738(2) 


68.513(2) 


48.1 


8 


36.9335(3) 


36.6235(3) 


23.9 


9 


24.58761(8) 


24.20344(7) 


14.9 


10 


18.38543(2) 


17.93780(2) 


10.6 


15 


8.6540358(4) 


7.9964587(2) 


4.2 


20 


6.13215967(2) 


5.36076877(1) 


2.7 



Table 2: Comparison with the Monte Carlo simulations by Janke and Kappler(1997). 







q = 10 


g = 15 


g = 20 


Cd 


large-g 
Monte Carlo 


18.38543(2) 
18.44(4) 


8.6540358(4) 
8.651(6) 


6.13215967(2) 
6.133(4) 


Co 


large- g 
Monte Carlo 


17.93780(2) 
18.0(1) 


7.9964587(2) 
7.99(2) 


5.36076877(1) 
5.361(9) 



g = 5 where the correlation length is as large as 2500. As for the asymptotic 
behavior of F^") at q ^ 4+, the Pade data oi F^ ' jx and Fq /x have the errors 
of a few percent around q = 4 and their behaviors are enough to convince us 
that the conjecture (|l|) is true for n = 2 with a — 0.073 ± 0.002 . Furthermore 

from the conjecture (|l|) the combination {F (n — |) |F("'|/F (|) F^^^} "~^ x~2 
is expected to approach the constant B for each n(> 3), and in fact the Fade 
data for every n{— 3, • • • , 6) gives B = 0.38 ± 0.05 , which also gives strong 
support to the conjecture for n > 3. 



4 Magnetization cumulants 

The behavior of the n-th magnetization cumulants M^"^ for /? ^ /3j at g = 4 



is well known as M 



(") 



A^:i{^Y 



and parallel to the conjecture for the 



energy cumulants by Bhattacharya et al. we can make a conjecture that 



M 



(") 

d,o 



(n) Mr 



(2) 



in the limit g — s- 4+ with (3 = (3t- We have examined the Fade approximation 
of M'^'^^D' for the large-g series generated by the finite lattice method for 
n = 2 and 3 both in the ordered and disordered phases, which in fact gives 



Table 3: The magnetic susceptibility for some values of q. 



q 


Xd 


Xo 


5 


9.13(3) X 10^ 


9.01(3) X 10* 


6 


6.585(4) X 10^ 


6.665(4) X 10^ 


7 


70.54(1) 


77.31(1) 


8 


19.359(1) 


21.525(1) 


9 


8.0579(1) 


9.0106(2) 


10 


4.23276(2) 


4.73823(4) 


15 


0.7304214(1) 


0.8056969(2) 


20 


0.309365682(1) 


0.33556421(1) 



Table 4: Comparison with the Monte Carlo simulations by Janke and Kappler(1997). 







q = 10 


g = 15 


9 = 20 


Mi^) 


large-q 
M.C. 


4.23276(2) 
4.233(2) 


0.7304214(1) 
0.73039(8) 


0.309365682(1) 
0.30936(4) 


M^^> 


large-q 
M.C. 


4.73823(4) 
4.74(4) 


0.8056969(2) 
0.805(3) 


0.33556421(1) 
0.3355(8) 



quite convergent Pade approximants for p — 15n/4 — 4 and as p leaves from 
this value the convergence of the approximants becomes bad rapidly again. In 
Table 3 we present the resulting estimates of the magnetic susceptibility Xd o — 
M^ o • ^^^ result is much more precise than the Monte Carlo simulatiortj at 
least by a factor of 100 as in Table 4. From the behavior of M^/"^ /x'-^^"'/^-'^'^ 
we obtain the coefficients in Eq. 



,(2) 



as ^i^^ = 0.0020(2), fi'^o"' = 0.0016(1) 
and fi'f'^ = 7.4(5) x 10-^ ^i^^ = 7.9(2) x 10"^ These convince us that the 
conjecture made for the magnetization cumulants is also true. 



5 Exponential correlation length 

Next we investigate the exponential correlation length ^i_o in the ordered phase. 
Here the exponential correlation length is defined by ^i = log (Ai/Aq) with Aq 
and Ai the largest and the second largest eigenvalues of the transfer matrix, 
respectively. (We have omitted the subscript T' in the previous sections) There 
is an obstacle to extract the correction to the leading term of the large-g expan- 
sion for the correlation length at the phase transition point from the correlation 
function < O(t)O{0) >c, since we know from the graphical expansion that it 



behaves like 

<O{t)O{0)>cOiz\l + 2zt^ + ---) (3) 

7^ exp {—mt). 

for a large distance t. Thus we will diagonalize the transfer matrix T directly 
for large-g. The eigenfunction for the largest eigenvalue Aq in the leading order 
of z is 

|0>= I P .0 > 

N 

where all of the N spin variables (each of which can take the value of s ~ 
0, 1, • ■ • , g — 1) are fixed to be zero, with the element of the transfer matrix 
< 0|r|0 > = 1 + 0{z'^) and the corresponding eigenvalue is Ag = 1 + 0{z). 
The eigenfunctions for the second largest eigenvalue are 

|/>= VlO ...0 eee ...e eO ...0 > 

(/ = 1,...,A^), 

9-1 



--;^l.l» 



.7-1^ 
with the diagonal matrix elements 

</|T|/>=z + 0(z2), 

and the off-diagonal matrix elements starting from higher orders in z. The 
second largest eigenvalues of T degenerate in the leading order with A^ = 
z-\-0{z^/'^). The degeneracy of the eigenvalues of the first A^ 'excited states' is 
the reason why the expansion series (Bf) of the correlation function cannot be 
exponentiated into a single exponential term. The off-diagonal matrix elements 
resolve the degeneracy with 

Ai/Ao = z + 4z3/2 + 6z2 + 0(z3), (4) 

Aat/Aq ^ z- 4z3/2 -f 6z2 -f 0{z^). 

for N —f oo. These eigenvalues constitute a continuum spectrum. It appears 
to be kept in any higher order of z. From Eq.(|j) we obtain 

l/6,o == - log z - 4zi/2 + 2z + 8/3x^/2 _^ (9(^2). 

6 



This ia.the same as the large- g expansion of 1/S,i,d given by BufFenoir and 
Wallom to this order. 

In the disordered phase the situation is quite similar. The eigenvalues 
of the transfer matrix for the first TV excited states constitute a continuum 
spectrum with their values exactly the same as in the ordered phase at least 
to the order of z^/"^ . 

The eigenvalues form the continuum spectrum just on the first order phase 
transition point. Off the transition point, we can see that the spectrum is 
discrete with A^ - A^+i ~ 0(e) for ^ <C e < 1 and A^ - A.^+i -- 0{e^/^) for 
e < V^ < 1 where e = (3/l3t - 1. 

6 Second moment correlation length 

Here we give the results of the large-g expansion of the second moment corre- 
lation lengths ^2nd,o in the ordered phase and ^2nd,d in the disordered phase 
defined by 

where /i2 and /xq are the second moment of the correlation function anrLthe 
magnetic susceptibility, respectively. The obtained expansion coefRcicnt£ll for 
the ordered and disordered phases coincide with each other to order z^ and 
differ from each other in higher orders. The ratio of the second moment corre- 
lation length in the ordered phase to that in the disordered phase is estimated 
by the Pade analysis to be not far from unity even in the limit of q ^ 4 
{S.2nd,o/S.2nd,d = 0.930(3)). Another point is that the ratio S,2nd,d/£.i,d of the 
second moment correlation length to the exponential correlation length is much 
less than unity in the region of q where the correlation length is large enough. 
It approaches 0.51(2) for q ^ 4. It is known that in the limit of the large 
correlation length. 



ei YZlcm^lil) 



<i, 



with ^i = — log(Ai/Ao). If the 'higher excited states' {i — 2,3, •••) did not 
contribute so much, this ratio would be close to unity, as in the case of the 
Ising pwd-cl on the simple cubic lattice, where £,2nd/£,i — 0.970(5) at the critical 
pointE3'E3 Our result implies that the contribution of the 'higher excited states' 
is large in the disordered phase of the Potts model in two dimensions even 
when q is close to 4. This strongly suggests that the eigenvalues of the transfer 
matrix for the first N excited states in the disordered phase form the continuum 
spectrum not only in the large-g region but also when q approaches 4. 

7 



As for the exponential correlation length ^i^o in the ordered phase for 
q ^ A, it is difficult to calculate the eigenvalues of the transfer matrix in 
much higher orders. It is quite natural, however, to expect that the ratio 
£,1.0/^1. d would be close to unity even in the limit of g ^ 4. The reason is the 
following. If the ratio ^i,o/$i,d would be 1/2 in the limit of g ^ 4, which is the 
known ratio in the second order phase transition point (g < 4), then the ratio 
S,2nd.o/S,i.o should bc closc to unity, which would imply that the higher excited 
states would not contribute so much to S,2nd.o and the eigenvalue of the transfer 
matrix for the first excited state would be separated from the higher excited 
states. This scenario is not plausible, since as already mentioned in section 5 
the continuum spectrum of the eigenvalues of the transfer matrix appears to 
be maintained in any high order in z. In this case, we can expect that the 
ratio ^2nd, 0/^1,0 would be around 1/2 as is the case in the disordered phase, 
resulting that ^i,o/Ci,d is close to unity. 



7 Summary 

We generated the large-g series for the energy and magnetization cumulants at 
the first order phase transition point of the two-dimensional g-state Potts model 
in high orders using the finite lattice method. They gave very precise estimates 
of the cumulants for g > 4 and confirmed the correctness of the Bhattacharya 
et al.^s conjecture that the relation between the cumulants and the correlation 
length for g = 4 and (3 ^> f3t (the second order phase transition) is kept in their 
asymptotic behavior for g — > 4_(. at /3 = /3i (the first order transition point). If 
this kind of relation is satisfied as the asymptotic behavior for the quantities 
at the first order phase transition point in more general systems when the 
parameter of the system is varied to make the system close to the second order 
phase transition point, it would serve as a good guide in investigating the 
system. 

Further the large-g expansion of the eigenvalues of the transfer matrix was 
calculated in the first 4 terms. We found that they have the same spectra of 
the eigenvalues of the transfer matrix in the ordered and disordered phases 
giving the same exponential correlation length (^1,0 = £,i.d) to the order of 
z^''^ and that the spectra are continuous in the thermodynamic limit. We also 
calculated the large-g expansion of the second moment correlation length in 
the ordered and disordered phases in high orders and found that they differ 
from each other in higher orders than z'^, but that the ratio ^2nd.d/£,i,d is not 
far from unity for all region of g > 4. We also found that S,2nd,d/S,d.i is far from 
unity even in the limit of g ^ 4. It receives significant contributions not only 
from the 'first excited state' but also 'higher excited states' and this suggest 



strongly that the continuum spectrum would be maintained (i.e. there would 
be no particle state) in the disordered phase. From these results it is quite 
natural to expect that the exponential correlation length ^i^o in the ordered 
phase is not far from that in the disordered phase even in the limit of q -^ 4 
and it is not plausible that their ratio approaches 1/2 that is their ratio in the 
second order phase transition point (q < 4). 
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